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Introduction: Geochemical models of core for- 

mation are commonly based on core and mantle abun- 
dances of siderophile elements that partitioned be- 
tween silicate and metal in a magma ocean in the early 
Earth. Tungsten is a moderately siderophile element 
that may provide constraints on the pressure, tempera- 
ture, composition, and oxygen fugacity conditions, and 
on the timing of core formation in the Earth. Previous 
experimental studies suggest that pressure exerts little 
to no influence over W metal-silicate partitioning up 
to 24 GPa 1 ' 8 , and indicate that the stronger influences 
are temperature, composition, and oxygen fugacity. 
However, core formation models based in part on W, 
predict metal-silicate equilibration pressures outside 
the available experimental pressure range, requiring 
extrapolation of parameterized models 3,5,6,7 . Therefore, 
higher pressure experimental data on W were needed 
to constrain this important parameter. 

Experiments: We measured the partitioning be- 
havior of W in high pressure and temperature experi- 
ments using the multi-anvil press (MAP) in the high 
pressure lab at NAS A- JSC, and using laser-heated di- 
amond-anvil cells in the Laboratory for Mineral Phys- 
ics at the University of Maryland (UMD). Starting 
materials consisted of a mixture of natural peridotite 
and metal powders with approximately 2:1 and 5:1 
wt% silicate to metal ratios for MAP and LHDAC ex- 
periments, respectively. The metals consisted of Fe and 
W in a ratio of approximately 5:1 wt%. 

The MAP experiments were conducted using 
COMPRES standardized octahedral assemblies with 
five experiments in MgO capsules and three in graph- 
ite capsules. The experimental run temperatures ranged 
from approximately 2400 to 2700 K, and pressures 
ranged from 10 to 26 GPa. 

LHDAC experiments consisted of a mixed silicate 
and metal foil that was approximately 5-8 microns 
thick, sandwiched between two insulating peridotite 
layers. Pressure calibration was by ruby fluorescence. 
Temperatures were determined by fitting measured 
thermal emissions to the Planck radiation function. 
Two experiments were at 3800 K and a third at 4400 
K. Pressures ranged from 28 to 50 GPa. LHDAC sam- 
ples were recovered for compositional analysis using a 
FEI Quanta 3D FEG focused ion beam (FIB) at 
NASA-JSC, yielding sample sections approximately 
one micron thick (Figure 1). 
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Figure 1. LHDAC sample G36; P = 28 GPa and T = 3800 K; 
qs-quenched silicate melt, qm-quenched metallic liquid, fp- 
ferropericlase. 


Recovered samples were analyzed using electron 
microprobes at both NASA-JSC (Cameca SX-100) and 
UMD (JEOL JXA-8900R). 

Results and Discussion: Liquid metal/liquid sili- 
cate partition coefficients (D w ) were calculated for 
each experiment using the mass fractions of W in the 
metal and silicate phases. Experimental oxygen fugaci- 
ty was determined from the Fe contents of coexisting 
silicate and metallic melts, assuming ideal solutions. 

Using the present data combined with literature da- 
ta i,2,4,5,6,7,8^ t j ie me t a l_ silicate partitioning of W was 
parameterized using multiple linear regression to yield 
the regression equation log D w (±0.30) = a + b (AIW) 
+ c/ T+ dP/T+e ZfiXt + fX c + g X s ; the silicate 
compositional parameters (fX t ) are the mole fractions 
of Si0 2 , A1 2 0 3 , CaO and MgO, and the metal parame- 
ters X c and X s are the mole fractions of C and S. Re- 
gression coefficients and their 1 -sigma uncertainties 
are: a = 1.07 (0.61), b = -1.41 (0.12), 

c = -2669 (550), d = -83 (11), e [Si0 2 = 2.49 (0.80), 
Al 2 0 3 = 5.99 (1.63), CaO = -10.4 (1.21), 

MgO = -3.10 (0.79)], f= 1.78 (0.33), and g = -2.75 
(0.78). 

Using this regression equation, log D w was calcu- 
lated along a peridotite liquidus 9,10 with f0 2 = -2.3 
AIW, silicate compositional terms of X S i 02 = 0-38, 
Xai 203 = 0.02, X Ca o = 0.03, and X Mg0 = 0.51, and metal 
compositional terms of X c = 0.009 and X s = 0.03. The 
result is that W is increasingly lithophile with increas- 



in g pressure along a peridotite liquidus, i.e. in a deep- 
ening magma ocean. Isolating the pressure and tem- 
perature dependencies by calculation along isobars and 
isotherms, log D w varies by approximately 1.4 log 
units over 0-60 GPa, and by approximately 1.0 log unit 
over 2000-4000 K. Within uncertainties, these rela- 
tionships indicate that the pressure and temperature 
dependencies are of roughly equal importance across a 
P-T range applicable to probable conditions of a deep 
magma ocean. 

Applying this parameterization, and literature pa- 
rameterizations 11 for Ni and Co, to the problem of core 
formation in the Earth, a pressure-temperature solution 
of 50 GPa and 3600 K was determined along a 
peridotite liquidus for a single-stage, magma ocean 
core formation model. A total pressure-temperature 
solution space of 48-54 GPa and 3500-4200 K was 
determined (Figure 2). A continuous accretion model 
based on exponential planetary mass growth rate, line- 
ar mantle and core density growth rates, and homoge- 
neous accretion resulted in a pressure-temperature so- 
lution of 59 GPa and 3900 K that represents the final 
depth of the magma ocean along a peridotite liquidus. 
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Figure 2. Pressure-temperature solution space for equilibrium 
core formation model. Colored regions are ranges of P-T 
conditions where the predicted partition coefficients equal 
core-mantle partition coefficients; widths of bands represent 
uncertainty on core and mantle abundances. Black polygon is 
the total solution space for W, Ni and Co. Liquidus is a fit to 
average of a peridotite liquidus and a chondritic mantle 
liquidus 9,10 . Core-mantle partition coefficients based on sili- 
cate Earth and core compositional models 12 ’ 13 ’ 14 . 

The timing of core formation in the Earth can be 
constrained using the Hf-W geochronometer. The 
magnitude of the W isotopic anomaly in this system is 
directly related to the partitioning of W into the core. 
Using the W parameterization of this study, the Hf/W 


ratio is predicted to decrease substantially with increas- 
ing pressure and depth. A primitive mantle Hf/W = 20 
is predicted along a peridotite liquidus at approximate- 
ly 53 GPa and 3800 K (Figure 3). 



Figure 3. Evolution of Hf/W ratio during core formation, 
calculated along a peridotite liquidus, assuming 60% volume 
fraction mantle melting and 3 mol% S in the core. 


Conclusions: This research has extended the ex- 
perimental pressure and temperature ranges of W met- 
al-silicate partitioning data and thereby provides better 
constraint on model predictions that previously re- 
quired extrapolation from lower P-T data. We find that 
the temperature and pressure effects on W metal- 
silicate partitioning are essentially equal. Single-stage 
and continuous accretion models based on W, Ni and 
Co show pressure and temperature ranges indicative of 
equilibrium core formation in a deep magma ocean; a 
Hf/W ratio that is consistent with a primitive mantle 
value was predicted at similar pressure-temperature 
conditions. 
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